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We examine dark energy models in which a plrantom field cj> is rolling near a local minimum of 
its potential V {(j)). We require that {l/V){dV/d4>) < 1, but {l/V){d^V / dffP) can be large. Using 
techniques developed in the context of hilltop quintessence, we derive a general expression for w as a 
function of the scale factor, and as in the hilltop case, we find that the dynamics of the field depend 
on the value of (l/V) [d^V/dfjP) near the mimimum. Our general result gives a value for w that 
is within 1% of the true (numerically-derived) value for all of the particular cases examined. Our 
expression for w{a) reduces to the previously-derived phantom slow-roll result of Sen and Scherrer 
in the limit where the potential is flat, {l/V){dV /dcj)) <C 1. 



I. INTRODUCTION 

Considerable evidence [l], Q has accumulated suggest- 
ing that at least 70% of the energy density in the universe 
is in the form of an exotic, negative-pressure component, 
called dark energy. (See Ref. Q for a recent review). 

A parameter of considerable importance is the equa- 
tion of state (EoS) of the dark energy component, defined 
as the ratio of its pressure to its density: 
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(1) 



Observations constrain w to be very close to —1. If ui is 
assumed to be constant, then —1.1 ^ —0.9 On 
the other hand, a variety of models have been proposed 
in which w is time varying. A common approach is to use 
a scalar field as the dark energy component. The class 
of models in which the scalar field is canonical is dubbed 
quintessence d, 0, H, S [l3| and has been extensively 
studied. 

A related, yet somewhat different approach is phantom 
dark energy, i.e., a component for which w < — 1. The 
simplest way to realize such a component is to use a scalar 
field with a negative kinetic term in its Lagrangian, as 
first proposed by Caldwe ll llll . Such models have well- 
known problems [H, [13, fl^ T 15 1 , but nevertheless have 
been widely studied as potential dark energy candidates 

[11 [11, [II [11, [11, [13, El, [12, [11 [11 ii] . 

Given the considerable freedom that exists in choosing 
the potential function of the scalar field, it is useful to 
develop general expressions for the evolution of w which 
cover a wide range of models. One such approach is to be- 
gin with the observational result that w for dark energy is 
very close to —1 today, and to assume that w was always 
close to —1 in the redshift regime of interest. With this 
assumption, one can assume an expanding background 
that is very close to ACDM and solve for the evolution of 
the scalar field for this case. This assumption alone is not 
sufficient to derive a general solution for the evolution of 
the scalar field or its equation of state. However, certain 
fairly general classes of such models can be solved exactly. 
The simplest such models, explored in Ref. [11|, assume 
a scalar field initially at rest in a potential satisfying the 



"slow-roll" conditions: 
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The first condition insures that w is close to — 1, 
while the two conditions taken together indicate that 
(lIV Mdy I d4)) is nearly constant. In the terminology of 
Ref. [23|, these are "thawing" models. 

For all potentials satisfying these conditions, it was 
shown in [26| that the behavior of w can be accurately 
described by a unique expression depending only on the 
present-day values of 17^ and the initial value oi w. In 
[28| this result was extended to phantom models satisfy- 
ing Eqs. and the w dependence of these phantom 
models was shown to be described by the same expression 
as in the quintessence case. 

The slow roll conditions, Eqs. (HJSl), while sufficient 
to ensure w ~ — 1 today, are not necessary. In [1^, a 
second possibility was considered, in which equation ([2]) 
holds, but equation ([3]) is relaxed. This corresponds to 
a quintessence field rolling near a local maximum of its 
potential. As in the case of slow-roll quintessence, this 
case can be solved analytically. In this case, there is an 
extra degree of freedom, the value of {l/V)(d^V/ d4>^), so 
that instead of a single solution for the evolution of w, one 
obtains a family of solutions that depend on the present- 
day values of fi^ and w and the value of {1 /V){d'^V / d^P') 
at the maximum of the potential. This family of solutions 
includes the slow-roll solution as a special case in the limit 
where {l/V){d'^V/d(l)'^) 0. 

Here we complete this series of studies by extending 
the above result to phantom dark energy models. Since 
phantom fields roll up their potentials, the analogous sit- 
uation is a phantom rolling close to a local minimum in 
its potential. We find that a unique family of solutions, 
very similar to the one derived in Ref. [1^, can be used 
to approximate the behavior of w in these models. 
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II. PHANTOM EVOLUTION NEAR A 
MINIMUM 

In this section we derive a general expression for the 
evolution of w for a phantom field near its minimum. Our 
treatment is similar to [1^. First consider a minimally 
coupled phantom field in a potential V {4>) . The phan- 
tom field has a negative kinetic term in its Lagrangian, 
leading to an equation of motion 



Substituting the above expansion into Eq. Q , and taking 
V {(j)*) — p^o we obtain 



il - [V" + {3/A)V {(l),)]u = 0. 
With the definition 



(4) 



fc = + (3/4) 1/(0,), 

we obtain the general solution to Eq. (fTTj) to be 
u = A sinh {kt) + B cosh (kt) . 



(11) 
(12) 

(13) 



where a is the scale factor and H = d/a is the expansion 
rate. Dots denote derivatives with respect to time and 
primes denote derivatives with respect to the field 4>. In 
a flat universe, the expansion rate is linked to the total 
density px via the Friedman equation (in units where 
8ttG = 1) as 



= Pt/3. 



The evolution of the scale factor is given by: 

Pt) , 



a 
a 



-lipT 



where is the total pressure 
Using the transformation 



0(0 = w(^)/a(^)^/^ 



Eq. (HI becomes 
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Now consider a universe consisting of pressureless mat- 
ter and a phantom field, where the phantom plays the 
role of the dark energy. In order to realistically mimic 
the observed dark energy, the phantom must have w close 
to —1 and its energy density must be roughly constant. 
The total pressure is then simply given by px ~ —PipOj 
where p^o is the present day density of the dark energy. 
(In what follows, a subscript always indicates a present 
day value). Under this approximation, Eq. ([H]) becomes: 



(9) 



We now apply Eq. ([9]) to a phantom rolling near a local 
minimum 0, in its potential. For any close to the min- 
imum, the potential can be expanded in a Taylor series: 

F (0) = 1/ (0*)+(l/2) V" (0,) (0 - 0O'+O ((0 - 00 

(10) 



The requirement of w « — 1 implies that the potential 
term dominates the kinetic term in the phantom's energy 
density. In such a scenario the evolution of the scale 
factor can be well-approximated by ACDM (see e.g. [soj): 
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sinh2/3 (t/tA) , (14) 



where fi^o is the present-day value of the phantom den- 
sity parameter fi^ and a = 1 at present. The time t\ is 
defined as 



tA = 2/^/3^ = 2/V3F(0,). 



(15) 



These results give the general solution to equation 
Eq. (H]) (under the approximations described above) as: 



0(0 = 



1 - n 



4>o 



1/2 



A sinh (kt) + B cosh (kt) 
sinh(t/tA) 



(16) 



where A and B are arbitrary constants. If we require 
(t = 0) = 0j, then B = and 



0i sinh {kt) 
ktA sinh {t/tA) 



(17) 



This is identical to the evolution of 0(0 for quintessence 
near a local maximum in the potential [29| . 

The EoS parameter w for a phantom is given by 



1 + w 



(18) 



Our requirement that w ~ —1 implies that w p^Q w 
F(0*). This, together with Eq. ^ and Eq. ^ gives 
an expression for w{a) (normalized to wq, the present-day 
value of w): 
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1 + w{a) = (1 + Wo)a 



/rj^ofc^A cosh [kt (a)] — ^(1 — f20o)a ^ + ^it>o sinh (a)] 



/n^o^^A cosh(fcio) ^ sinh(fcio)] 
I 



(19) 



where t{a) and can be derived from Eq. H 



and 



t{a) ~ sinh 



to — t\ tanh ""^ 



1 - r2 



00 



(20) 



(21) 



For convenience, we now switch to the constant K = 
ktA. In terms of the phantom potential, K can be written 
as 



J 



= v/i + (4/3)y"(0,)/y(0,). 



(22) 



indicating that K (which is always > 1 since V (0* ) > 
and V" (0*) > 0) depends only on the value of the 
potential and its second derivative at its minimum. In 
terms of K we can express the evolution of w in the 
following form: 



, _ , ^ ^„3(K-l) ma) + 1)^ {K F{a)) + (f (a) - !)'< {K + f (g))]^ 



r 



(23) 



where F{a) is given by 



F{a) = ^l + {n-^-l)a- 



(24) 



(Note that F{a) — 1/ y/n^(a), where fl^{a) is the value of 
as a function of redshift, so that F{a = 1) = fl^^ .) 

Eq. is our central result. We therefore find that, 
in the limit where w is slightly less than —1 (i.e., where 
the phantom potential energy dominates its kinetic en- 
ergy), all phantom models with a given value of wq and 
n^o converge to a unique evolution of w characterised 
by the value of K (which depends on V" {(j)^) /V^(0,), 
i.e., loosely speaking, on the curvature of the potential 
at its minimum) . For phantom potentials which are very 
flat, i.e., which satisfy both slow-roll conditions, a simi- 
lar analysis was shown to yield only a single form of w{a) 
(Eqn. (18) of nil). Eq. (|^ captures the more complex 
behavior introduced by the dependence of w{a) on K. It 
is straightforward to show analytically that Eq. (|23p goes 
over to Eqn. (18) of [1^ in the slow-roll limit, K ^ 1. 

Note that Eq. (|23|) has an identical functional form to 
the corresponding expression for w{a) for quintessence 
models derived in Ref. [2^. The only difference for 
the phantom models is that 1 + wo < 0, producing a 
corresponding w{a) < —1. Of course, this is a signifi- 
cant difference in deriving observational constraints on 
the models. Further, there is a sign difference between 
the definition of K for phantom models (Eq. and 



for quintessence models [29j, corresponding to the fact 
that here we are dealing with a minimum in the pote- 
nial, rather than a maximum. 

Interestingly, Eq. reduces to simple polynomial 
forms if if is a small integer. Some of these forms are 
given in 29]. While it would require fine tuning of the 
potential for K to be equal to a small integer, these forms 
are still useful in obtaining qualitative insight into the 
behavior of w (a) as a function of K. 



III. COMPARISON TO EXACT SOLUTIONS 

We now turn to comparing our analytic expression for 
the evolution of w to the numerically computed exact 
evolution for a few different models. In each case we 
have a perfect fluid dark matter and a phantom fleld 4> 
dynamical dark energy. 

We now consider three different phantom potentials 
which have local minima. The phantom analog of the 
PNGB model [3l| , has a potential given by 



n0)-P0o + M^ [l-cos(0//)] 



(25) 



where M and / are constants. Other models with a local 
minimum in the potential include the Gaussian potential. 



F(0) = p^o + 



1 



(26) 
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and the quadratic potential 



(27) 



where a and V2 are constants. 

We set initial conditions deep within the matter- 
dominated regime. We choose wq = —1.1 so that our 
results will give an upper bound on the error in our ap- 
proximation for —1 > Wq > —1.1. The value of the 
potential at the minimum, y (0*), is chosen to be equal 
to the energy of the cosmological constant. The initial 
value of the field (pi is taken to be slightly displaced from 
its minimum and is fixed to give wq — —1.1. The 
initial velocity of the field is taken to be zero. 

As discussed above, our formalism applies to phantom 
models for which Eq. ([2]) is satisfied, but Eq. ^ is not. 
The latter clearly holds when iiT ~ 1 but is violated as 
K departs from 1. For our specific examples, we focus 
on the cases K = 3 and K = 4. For K = 3, for all the 
potentials, 



1^' 

1 (PV 



O [10-1] 



V 



a— >0 



a— ►0 



0[1] 



For K = A, for all the potentials, 



1 dV 
1 d^V 



V 



a— ►0 



a— »0 



O [10-2] 
O[10] 



In Figs. ([T][3]), the evolution of w from Eq. P5)) is 
shown in comparison to the exact evolution for the three 
different models. The agreement, in all three cases, be- 
tween Eq. psp and the exact numerical evolution is ex- 
cellent, with errors Sw < 0.2% for the quadratic poten- 
tial, Sw < 0.6% for the PNGB potential, and 5w < 0.8% 
for the Gaussian potential. 

We note that the accuracy is larger for the quadratic 
potential than it is for the other two potentials. This 
is expected, since the derivation of Eq. (|23|) was based 
on the quadratic potential. This is however not the case 
for hilltop quintessence (29j . where in spite of a simi- 
lar derivation, the highest accuracy was for the PNGB 
potential. This is most likely due to an accidental can- 
cellation of errors. 

Finally, in Figs. [4] and \E\ we use our approximation 
(Eq. to construct a likehhood plot for wq and 

il^o with K — using the recent Type la Supernovae 
standard candle data (ESSENCE-f SNLS+HST from 0]). 
Clearly, these models are not ruled out by current super- 
nova data. As in the case of hilltop quintessence [29j . 
we find that a larger K increases the size of the allowed 
region. 

Eq. therefore allows one to map the w{a) behav- 
ior of a wide range of phantom dark energy models on 




FIG. 1: Comparison between our approximation for w{a) 
(Eq. I23|) with wo — —1.1 and il^o = 0.7 and the exact 
(numerically-integrated) evolution for w{a) for the quadratic 
potential (Eq. (|27|l '). Red (dashed) curves give our approxi- 
mation, and solid (blue) curves give exact evolution, for (left 
to right), a: = 3, 4, where K is defined by Eq. 1(22^ . 




FIG. 2: As Fig. 1, for the PNGB potential (Eq. (l25ll l. 




FIG. 3: As Fig. 1, for the Gaussian potential (Eq. (gSJ). 
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FIG. 4: Likelihood plot from SNIa data for the parameters 
Wo and 0,4,0, for phantom models with generic behavior de- 
scribed by Eq. (|23|l . with K = 3, where K is the function of 
the curvature of the potential at its maximum given in Eq. 
(|22[) . The yellow (light) region is excluded at the 2a level, and 
the darker (orange) region is excluded at the la level. Red 
(darkest) region is not excluded at either confidence level. 




FIG. 5: As Fig. 4, for K = 4. 



parameterization of these models in terms of fl^o, wq and 
K. Such a parameterization is particularly interesting in 
the light of recent work implying that the CPL parame- 
terization does not satisfactorily fit the SN+BAO+CMB 
data simultaneously at both low and high redshifts [34] . 

Subsequent investigations have shown that Eq. ([^ 
has a much wider applicability than the models described 
in this work and [29|. It has been found to describe 
quintessence (phantom) models in which the field rolls 
into a minumum (maximum) of the potential, leading to 
a richer set of behaviors including oscillatory solutions 
[S^ . In [3^ it was shown that Eq. I|23p is not limited 
to models evolving near extrema of their potentials and 
can be applied to a broader class of slow-roll quintessence 
models. [37] showed that Eq. ([23]) also applies to some 
slow-roll k-essence models. 



IV. CONCLUSION 

Using techniques previously applied to quintessence, 
we have derived a general expression for the evolution 
of w, which is valid for a wide class of phantom dark 
energy models in which is the field is rolling close to a 
local minimum. Such models provide a mechanism to 
produce a value of w that is slightly less than —1. We 
have tested our expression against the (numerically de- 
termined) exact evolution for three different models and 
in each case it replicated the exact evolution studied with 
an accuracy greater than 1%. A comparison between our 
generic approximation and the observational data indi- 
cates that these models are allowed by SNIa data, and 
that the size of the allowed region increases with the cur- 
vature (V" /V) of the potential at the minimum. 



to a common model-independent evolution. Clearly, this 
result is not well described by popular linear (CPL) pa- 
rameterization of w ^32., ,33], i.e., w{a) — wq + 'Wi{l — a), 
except in the special case of if — > 1. However, since our 
treatment demonstrates that the evolution of w is com- 
pletely generic to any model satisfying Eq. and de- 
pends only on fl^o, wq and the curvature of the potential 
at its minimimum (K), Eq. (|23p allows for an excellent 
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